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The Einstein-Cartan-Sciama-Kibble theory of gravity removes the constraint of general relativity 
that the affine connection be symmetric by regarding its antisymmetric part, the torsion tensor, as 
a dynamical variable. The minimal coupling between the torsion tensor and Dirac spinors generates 
a spin-spin interaction which is significant in fermionic matter at extremely high densities. We show 
that such an interaction averts the unphysical big-bang singularity, replacing it with a cusp-like 
bounce at a finite minimum scale factor, before which the Universe was contracting. This scenario 
also explains why the present Universe at largest scales appears spatially flat, homogeneous and 
isotropic. 
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The Einstein-Cartan-Sciama-Kibble (ECSK) theory of gravity, like general relativity (GR), is based on the gravita- 
tional Lagrangian density that is proportional to the curvature scalar R [I]. It removes, however, the GR constraint 
that the affine connection 1^* be symmetric by regarding the antisymmetric part of the connection, the torsion ten- 
sor S k i j — I^*], as a dynamical variable [2]. Varying the total Lagrangian density —^R\f—g + £m> where £ m is 
the Lagrangian density of matter, with respect to the contortion tensor CV^ = <Sijfc + Sjfci + Skji gives the Cartan 
equations 

S\ k -S l 5i + S k 5l = - ] jKSl ^ (1) 

where S% = S ik and s l ^ k = 2(5£, m /5Cijk)/ \J—g is the spin tensor. These equations are linear and algebraical: torsion 
is proportional to spin density and vanishes outside material bodies. 

Varying the total Lagrangian density with respect to the metric tensor gi k gives the Einstein equations with terms 
on the curvature side that are quadratic in the torsion tensor. Substituting ([l]) into these equations leads to the 
Einstein-Cartan equations dk — n(Tik + Uik), where dk is the Einstein tensor, T ik = 2(5£ m /5g th )/ V~5 is the 
energy-momentum tensor, and 

TT ik — k( — « ij s kl — -Ji l v k _|_ ioJ'K, k _|_ - n ik ( — A? 1 e JTO -l_«J 7m «., n P)\ 
U — K\ 8 [; S -S S 3l + -S S jl + -g ( 4S j[m S ;] +S Sjlm) J ( z ) 

is the correction to the energy-momentum tensor generated by torsion and quadratic in the spin tensor [3] . The spin 
tensor also appears in Tjfc because £ m depends on torsion. The corrections from the spin tensor to the right-hand 
side of the Einstein equations are significant only at extremely high densities, on the order of the Cartan density [4]. 
Below this density, the predictions of the ECSK theory do not differ from the predictions of GR. In vacuum, where 
torsion vanishes, this theory reduces to GR. The ECSK gravity therefore passes all observational and experimental 
tests of GR 0. 

In the ECSK theory, the Dirac Lagrangian density for a free spinor if) with mass m, minimally coupled to the 
gravitational field, is given by £ m = i,\/—g(ipY4';i ~ , 4 , ;iJ l '4 > ) ~ 9'4 ! '4 ) , where g is the determinant of the metric 
tensor and 7* are the Dirac matrices obeying = g lk I. Semicolon denotes a covariant derivative with respect to 

the affine connection [21 [3] : 

Ak = ^:k + \c ij kl [l l 3 H, = i>-.k- jCijfc^y 1 . (3) 
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and colon denotes a Riemannian covariant derivative with respect to the Christoffel symbols. We use the units in 
which c = fi = fcB = l, sok = 
is completely antisymmetric: 



which c = h = fee = 1, so k = SttG = m p 2 , where mp is the reduced Planck mass. For a Dirac field, the spin tensor 



s ijk = _ e m su a^ifyVV, (4) 

where e % i kl — e lJ /y/—g, e ij is the Levi-Civita permutation symbol, and s l is the Dirac spin pseudovector. The 
Cartan equations for such a field give therefore the completely antisymmetric contortion tensor [2J: 

Cijk — Sijk = —KCijklS ■ (5) 



Substituting Q into Q gives 



U lk = ^n(2s l s k + s l s ig lk ). (6) 



Varying £ m with respect to the spinor adjoint conjugate -0 gives the Dirac equation i^tp-^k = mip, which can be 
written as vy i/j-k = mip — |'«(?/'7 fc 7 5 '0)7fc7 5 ?A and whose conjugate is — i"4>;kl k — mrf ~ § K (?/'7 fe 7 5 ' ! /')V'7fc7 5 [31 H3- 
Putting this equation in the energy-momentum tensor corresponding to the Dirac Lagrangian, Tik — §(V ,< ^(i7fc)V';;; ~ 
$\j 6 ( t lk)^) ~ - 4>;jl J ^)9ik + mijipgik, gives 



1 _ _ . 

T ik = £ (#(,7*)^ - ^;3 S (ilk)^)- (7) 



Substituting ([3| and §5§ into ^ gives 



j _ . _ 1 

T ik = 2"(^ 5 (i7fe)^y - ip-.j^ff^ip) + -n(-SiS k + s l sig ik ). (8) 

The combined energy-momentum tensor for a Dirac field on the right-hand side of the Einstein-Cartan equations is 
thus 

i - ■ - 3 

Tik + Uik = A^lftk^-.j ~ ^:3 S [i1k)^) + ^s l sig ik , (9) 

which has been found by Kerlick [B]. This expression agrees with 0, where we derived it from the Hehl-Datta 
equation, which is the Dirac equation upon substituting Q and ^ 5J. The first term on the right of ^ is the 
GR part of the energy-momentum tensor for a Dirac field and can be macroscopically averaged at cosmological scales 
as a perfect fluid with the energy density e and pressure p. In the comoving frame of reference, in which go a = (a 
denotes space indices) and the four- velocity u l of the cosmological fluid satisfies u° = 1 and u a = [5J, the relation 
s l u.i = gives s° = 0. Accordingly, the second term on the right of (|9j) is equal to ~jns 2 gik, where s is the spatial 
spin pseudovector which measures the spatial density of spin. The average value of its square is (s 2 ) = |n 2 , where 
n is the fermion number density. The averaged second term on the right of ^ acts thus like a perfect fluid with a 
negative energy density: 

9 

e = — p = —an 2 , a = — k. (10) 
y ' 16 K ' 

Hchl, von der Heyde, and Kerlick have used the spin-fluid approximation of fcrmionic matter, Sijk — SijUk and 
SijU^ = 0, to show that the spin-density contribution to Tik + Uik behaves like a stiff matter with e = p = —jus 2 , 
where s 2 — \sikS lk = |n 2 [HI E| - This behavior is significant in spin fluids at extremely high densities, even without 
spin polarization, leading to gravitational repulsion and avoidance of curvature singularities by violating the energy 
condition of the singularity theorems [5]. Trautman, Kuchowicz, and others have shown that such a repulsion replaces 
the big-bang singularity with a nonsingular big bounce, before which the Universe was contracting |10LI11) . In contrast 
to spin fluids, Dirac spinors coupled to torsion enhance the energy condition for the formation of singularities [S"l 1 1 2 j . 

The spin-fluid model can be derived as the particle approximation of multiple expansion of the integrated con- 
servation laws in the ECSK gravity [T3]. The particle approximation for Dirac fields, however, is not self-consistent 
0]. The spin- fluid description also violates the cosmological principle [13]. In this paper, we use the Dirac form of 
the spin tensor for fermionic matter, Sijk — sum [2j |6j[T2], which follows directly from the Dirac Lagrangian and is 
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consistent with the cosmological principle [T5]. We show that the minimal coupling between the torsion tensor and 
Dirac fermions, despite enhancing the energy condition, also averts the big-bang singularity. 

As in |16j . we consider a closed, homogeneous and isotropic universe, described by the Friedman-Lemaitre- 



Robertson- Walker (FLRW) metric. In the isotropic spherical coordinates, this metric is given by ds = dt 



kr 2 /4) _2 (dr 2 +r 2 dd 2 + r 2 sin 2 1} dtp 2 ), where a(t) is the scale factor and k = 1. The corresponding Einstein equations for 
the combined energy-momentum tensor ([9| in the comoving frame become the Friedman equations (the cosmological 
constant is negligible in the early Universe): 



k = 



a 



3 

2ad + k 



-k(p 



+ an a 



(11) 
(12) 



where dot denotes differentiation with respect to the cosmic time t. These equations yield the conservation law 
Jj((e - an 2 )a 3 ) + (p + are 2 )|(a 3 ) = 0, which gives 



(13) 



and 



a A de - 2aa 6 ndn + (e + p)d{a A ) = 0. 

As in |17j . we use e, p and n for ultrarelativistic matter in kinetic equilibrium: e(T) = ^(/.(T^T 4 , p(T) 
n (T) = ^)-g n (T)T 3 , where T is the temperature of the early Universe [IS]. The effective numbers of thermal degrees 
of freedom are g*(T) — gb(T) + gfff(T) and g n (T) = jgf(T) (only fermions contribute to torsion), where g^ = YliQi 
is summed over relativistic bosons, gf — g t is summed over relativistic fermions, and is the number of the spin 
states for each particle species i. Substituting these values to (13) gives 



dT 



da 



= 0, 



(14) 



where h+ 
gives 



^(/^(T) and h n = t ^§rg n {T) can be assumed constant in the range of T considered. 1 Integrating (14) 



a r T r ( 3a ft. 2 



T 



-exp 



T~ 



(15) 



where a r is the scale factor at a reference temperature T r . 

The function a{T) (15) is not monotonic. As T increases, a decreases until T reaches a critical temperature T C] 
given by ^(T cr ) = 0, 



2ft, 
3aft 2 



1/2 



(16) 



and then increases. Since an increasing function a(T) is unphysical, a c 
the scale factor: 



= a t T T 



2ft. 



1/2 



a(T cr ) > is the smallest allowed value of 



(17) 



For T <§C T cr , (15) reduces to a — 93 ^ !L , which is satisfied in the 



The Universe is therefore nonsingular: a > a c 
radiation-dominated era. 

To verify that a cr is the minimum scale factor of the Universe, we substitute (14) into (11) without the negligible 
term k = 1, obtaining 



T 



1 



3aft 2 



T 2 2ft, 



'{KT 2 -ah 2 n T 4 



(18) 



1 For constant values of g t and g n , the relations e <x T 4 , p = | and n <x T 3 are consistent with a relation ^ = used in |16| 
e = p oc n 2 as in [l[9l[l6], then we also have ^ = ^j^pg ■ 
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Denoting /3 = T _1 and /3 cr = T" 1 leads, with rfisb , to 



3 Her 



/? 2 - /3 C 2 



(19) 



which yields /3 > /3 cr and T < T cr . Equation (151 gives then a > a cr . We make a substitution: 



/3 cr coshr/, 



(20) 



where 77 is a parameter satisfying \rj\ > r) CT = arcoshw ^. Putting (201 in (19) gives to(4cosh V ~ 1 



lfl = 



1, where 



to = /3 2 r J is a characteristic time scale of the torsion-dominated era. Integrating this equation, with t = taken 
as the instant at which 77 jumps from — 77 cr to ?7 C1 ., gives 



t 1 • , /„ N 2 

— = -smh 2?7 - -77 + 
to 6 3 

t 1 ■ , s 2 

— = -smh 2?/ - -77 - 
to 6 3 



v/3 2 

~eT ~ 3 



?7cr, ?7 < -T^cr, 



/3 2 
3 



(21) 



The parametric Eqs. (20) and (21) determine Pit). Putting f3(t) in (15), written by means of (17) as 



a r /3 



/3/ eXP V2^ 



Pi 



a C rP 



^/3 CT eXP V2/3 2 



B 2 

r'c 



(22) 



gives the dynamics of the early Universe, a(t). 

Figs. [T] and [2] show how the temperature and the scale factor, respectively, depend on the cosmic time. The 
bouncing point, where the scale factor has its minimal value a CT , is a cusp. As 77 increases from some initial, negative 
value, the Universe contracts (d < 0) until 77 = — ?7 cr , at which /3 = j3 CI and a — a cr . Then it undergoes a cusp-like 
bounce, at which 77 jumps from — 77 cr to ?7 cr . After the bounce, 77 increases to infinity and the Universe expands (d > 0). 
The unphysical big-bang singularity appearing in general-relativistic cosmology is replaced in the ECSK gravity by a 
nonsingul ar ( with respect to curvature) big bounce that follows a contracting phase of the Universe [11]. For /? 3> /3 cr , 
(19) and (22) give a oc T~ x oc t 1 / 2 , which is characteristic to the radiation-dominated era. 




FIG. 1: The normalized temperature 3- as a function of the normalized cosmic time 



The co ndit ion d = 0, defining the scale factor at a stationary state a = a s t, would be satisfied at a temperature T st 
jiven by ( 11 ) without the negligible term k = 1: /i*T 4 — ah 2 n T & = 0. This temperature is equal to 



A. 



1/2 



> T cr , 



(23) 
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FIG. 2: The normalized scale factor — as a function of the normalized cosmic time A. 



so the Universe never reaches T — T st and a = 0. 
bounce from a = — v to a = v, where 



At the minimum scale factor a cr , the Universe undergoes a cusp-like 



\a(T CI )\ 



QLfh~. J- „ r ) Ct „ , 



1/2 



32e 
243 



1/2 



K 
h„ 



a r T r . 



(24) 



At the bounce, the velocity of the point that is antipodal to the coordinate origin in a closed Universe is equal to 
fant(7ci') = ttv [If)] . The density parameter at the bounce is given by [TB] 



npy- 1 + ^-1+ 2uhl 



32e/i2(a r T r ) 2 

As the Universe expands, the antipodal velocity decreases according to 

1/2 

3 v— — „- - j " \ 

and the density parameter increases according to 



<\„niT) = rr./iri = ( ';'://. V' nh'iV'ur) ' rr( '>. /- n//;7 ; ;) ' a r T r exp^ 2 



n(T) = i 



i 



= 1 + (a r T r ; 



: (KT Z - ahiT^) ) exp - 



y2 
y2 



(25) 



(26) 



(27) 



(26) reduces to v ant (T) ~ T and (27 1 reduces to 



When the Universe reaches the radiation-dominated era, T -C Tor, 
O(T) - 1 - T- 2 . 

Kerlick has shown that Dirac fields in the ECSK gravity satisfy the strong energy condition for the Hawking- 
Penrose singularity theorems |19j . and that the torsion-induced spin-spin interaction (the second term on the right 
of (|9|) enhances that condition [B]. Such an interaction is thus attractive. O'Connell has obtained the same result 
using the spin contributions to the Lagrangian density [12] . Writing er^ = Tik + Uik, the strong energy condition 
holds if W = [oik — ^gik&\)'u 1, u k > 0. For a homogeneous Universe, in the comoving frame of reference of an 
irrotational cosmological fluid, the tensor (|9| and the Dirac equation give W = \m^}i\) J r |ks 2 , suggesting that torsion 
enhances the formation of singularities [6j . This result agrees with the contribution to W from the torsion-induced 
spin-spin interaction, e + 3p — 2an 2 > 0. The singularity theorems, however, have another important component: 
the Raychaudhuri equation, which describes the time evolution of the expansion scalar for timelike congruences |20| . 
The expansion scalar 8 — u l ^ measures the fractional rate at which a small volume of matter changes with time as 
measured by a comoving observer. Because of the Raychaudhuri equation, timelike geodesies for the matter satisfying 
the strong energy condition converge (8 decreases), resulting in a caustic, 8 — > — oo (if 8 decreases continuously), and 
thus in a singularity within a finite proper time |20j . 

For the cosmological fluid in the comoving frame, the FLRW metric gives 8 = (ln-y/— g) y o = • Since the Universe 
at the minimum scale factor a cr undergoes a bounce from a — — v to a — v, the expansion scalar has a discontinuity 
there, increasing from 8 ~ — — to 8 = — . As the Universe before the big bounce contracts, 8 decreases from some 



() 



value given by the initial condition of the contraction to a value — At a = a cr , 9 jumps to — . As the Universe 
after the bounce expands, 9 decreases towards zero. The discontinuity of the expansion scalar at the cusp-like big 
bounce therefore prevents 9 from decreasing to — oo and reaching a curvature singularity and guarantees that timelike 
geodesies in the Universe with Dirac fields coupled to torsion continue through the bounce. A similar discontinuity 
completes timelike geodesies in the gravitational field of an Einstein-Rosen bridge [2"T] . 

If we assume that the early Universe contained only known standard- model particles, then <?b = 28 and gt = 90 
[IB] , so the temperature at the big bounce (16) was 



T cr « 0.78 mp. 



(28) 



As the reference values, we can take the temperature and scale factor at the matter-radiation equality, where the 

6 cq 



radiation-dominated era ends. We therefore have T r = T cq ps 0.75 eV and a r = a oq = ^2" > where qq ps 2.9 x 10 27 m 



is the present scale factor [IB] and z cq ps 3200 is the redshift at the matter-radiation equality 18J. Accordingly, the 
values at the big bounce of the scale factor (17), antipodal velocity and density parameter were, respectively: 

-4 m 



,34 



a cr w 5.9 x 10 
«a„t(T cr ) ps 8.9 x 10 J . 
ft(T cr ) ss 1 + 1.3 x 10~ 70 . 



(29) 
(30) 
(31) 



The value ( 30 ) is enormous, whereas ( |31[ ) differs from 1 by a quantity that is extremely small in magnitude, like in 
the spin- fluid cosmology with torsion |16) . Such extreme values result from a r T r 3> 1. 

The big-bounce value of A = fl — 1 is on the order of | Q , where Qs is the spinor-torsion density parameter (16) . 
subscript R denotes radiation, and denotes the present value. An enormous value of a r T r is related to an extremely 
small magnitude of Qs'- 



n s A 



n 2 



(Gs 2 / Pc )A 



Gn 2 p c A 



GT 6 Pc A 



J>8 



H 2 A 



J>2 



(aTY 



[a r T r 



(32) 



where p c is the critical density and H is the Hubble parameter. The apparent fine tuning of Q(T CI ) is thus caused by 
|f^s| 1, originating from an extremely weak spinor-torsion coupling in the ECKS gravity, as in 16j. Accordingly, 
such a coupling naturally explains why the present Universe at largest scales appears nearly fiat, solving the flatness 
problem without introducing exotic matter fields necessary for cosmological inflation. 

This coupling is also responsible for an extremely rapid expansion of the Universe after the big bounce. Such an 



expansion produces an enormous number of causally disconnected volumes, iV ■ 



t , from a single causally connected 



region (the closed Universe before and at the big bounce), naturally explaining why the present Universe at largest 
scales appears homogeneous and isotropic. The spinor-torsion coupling therefore solves the horizon problem without 
inflation. The transition from the torsion-dominated era to the ra dia tion-dominated era occurs naturally as the 
contribution from this coupling to the Friedman equations (11) and (12) rapidly weakens (according to T~ 6 ), which 
is another advantage of this scenario. 
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